Helical edge and surface states in HgTe quantum wells and bulk insulators 
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The quantum spin Hall (QSH) effect is the property of a new state of matter which preserves 
time-reversal, has an energy gap in the bulk, but has topologically robust gapless states at the 
edge. Recently, it has been shown that HgTe quantum wells realize this novel effect In this 
work, we start from realistic tight-binding models and demonstrate the existence of the helical edge 
states in HgTe quantum wells and calculate their physical properties. We also show that 3d HgTe 
is a topological insulator under uniaxial strain, and show that the surface states are described by 
single-component massless relativistic Dirac fermions in 2 -|- 1 dimensions. Experimental predictions 
are made based on the quantitative results obtained from realistic calculations. 
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Conventional insulators have a gap for all charge ex- 
citations and their physical properties are not sensitive 
to changes in the boundary conditions. Recently, a new 
class of quantum spin Hall(QSH) insulators has been pro- 
posed theoretically in both 2d and 3dP, i, H i, Q. The 
QSH insulators are invariant under time reversal, have a 
charge excitation gap in the bulk, but have topologically 
protected gapless edge states that lie in the bulk insu- 
lating gap. This type of insulator is typically realized in 
spin-orbit coupled systems; and the corresponding edge 
states have a distinct helical property: states with one 
spin-polarization move around the sample edge in one di- 
rection, while states with the opposite spin-polarization 
move in the opposite direction [3, 0, 01- The helical edge 
states are responsible for the intrinsic spin Hall effect in 
the insulating stated, S HO, [HI- The edge states come 
in Kramers's doublets, and time reversal symmetry en- 
sures the crossing of their energy levels at special points 
in the Brillouin zone. Because of this level crossing, the 
spectrum of a QSH insulator cannot be adiabatically de- 
formed into that of a topologically trivial insulator with- 
out helical edge states; therefore, in this precise sense, 
the QSH insulators represent a topologically distinct new 
state of matter. In 3d the class of strong topological in- 
sulators has a similar distinction and is the natural 
generalization of the QSH insulator. 

The study of the QSH effect in quasi-2d HgTe/CdTe 
quantum wells carried out in Ref.jl'] is based on a sim- 
plified model obtained by k • P perturbation theory and 
the envelope function approximation. The conclusion of 



a topological quantum phase transition is reached based 
on a Dirac- type subband level crossing at the F point. 
However, such a level crossing is not generic and could be 
avoided in the real system due to the bulk inversion asym- 
metry (BIA) of the zinc-blende lattice. Consequently, a 
more realistic study is necessary to obtain a better un- 
derstanding of the QSH phase and the topological phase 
transition. In this paper, we study the subband struc- 
ture and edge state properties of HgTe/CdTe quantum 
wells using a realistic tight-binding(TB) model. The level 
crossing avoided at the F point is recovered at several fi- 
nite wavevectors. In other words, the phase transition 
between two different insulating regions remains robust 
despite inversion symmetry breaking. Furthermore, the 
topological nature of the QSH regime is demonstrated 
explicitly by studying the properties of the helical edge 
states in an open boundary system. We also apply the 
same realistic TB calculations to uniaxial strained 3d 
HgTe, and obtain the topologically non-trivial surface 
states. Thus, strained bulk HgTe is demonstrated to 
be a strong topological insulator. Ref. investigates 
a class of models with bulk inversion symmetry which 
strictly speaking does not apply to HgTe, however, adia- 
batic continuity is used to argue that strained bulk HgTe 
could be a strong topological insulator. 

The CdTe and HgTe materials have the same zinc- 
blende lattice structure and are well-described by the 
same type of tight binding Hamiltonian, albeit with dif- 
ferent parameters. This Hamiltonian includes two s-type 
orbitals and three p-type orbitals on each atom and reads 
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where and _> ^ are the creation operators for 

i(7 R 1(7, R-\- d 

electrons on the anion and cation sites respectively, i — 
{s, s* ,px,Py,Pz) is the orbital index, and a is the spin 
index. Ei^a, Ei^b and Vi_j are the tight binding parameters 
defined by Slater and Koster[l3|- Spin-orbit coupling is 
contained in the last two terms, and is represented by 
two coupling constants Ao and Af,. The tight binding 
parameters are taken from Ref. where they were 

determined by fitting to a first-principles calculation. 

To obtain the boundary (edge or surface) states for 
quasi-2d HgTe/CdTe quantum well structures and for 
the bulk 3d materials, we apply the Green's function 
method^, based on the TB model described above. 
For the quantum well system, we consider a symmetric 
HgTe/CdTe hetero-structure with a fixed Nc = 8 layers 
of CdTe surrounding a variable Nh layers of HgTe on each 
side. In order to calculate the boundary states we choose 
open boundary conditions along the a;-direction, periodic 
boundary conditions along the y-direction, and the open 
boundary condition along the z-direction which is the 
growth direction of the quantum well. For the bulk ma- 
terials, surface states can be calculated for a semi-infinite 
system with open boundary conditions in one direction, 
such as [001], and periodic boundary conditions in the 
other two directions. The inverse Green's functions of 
both the quantum well and bulk materials can be writ- 
ten in a block tri-diagonal form as 
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where the diagonal block describes the Hamiltonian 
within the same "principal layer", [l^ that is, the layer 
along the x-direction with open boundary conditions for 
the quantum well case, or that along the only direction 
with open boundary condition for the bulk case. The 
off-diagonal block C describes the coupling between two 
nearest-neighbor principal layers. To study the bound- 
ary states we only need g^, the Green's function on 
the boundary, where i, j are the indices of the local ba- 
sis on the boundary. This function is contained in the 
first diagonal block of the matrix G (z) = (z — 
and can be expressed in a recursive way as g. 
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{z-Ha- Cy^^-i^Ct) .^.^ with g'lp denoting the bound- 
ary Green's function for a system with TV principal lay- 
ers. The above recursive equations can be closed by 
the initial condition glj — (z — Hq) , and we ob- 
tain g^P"^ iteratively. Any physical observables pro- 
jected onto the boundary are easily expressed using these 
Green's functions as — ^ / du^^ - Img^^ {lo + iO~^)Oji = 

— i J dijjpo (uj) where po (w) is a type density of states and 
i,j run over all the local basis states on the boundary. 
For example, 0'= = Y.^ N >< and = J2ij Lfj + Sfj 
generate the density of states for charge and spin on the 
boundary respectively. Notice here we generalize the def- 



inition of "spin" to include all the local angular momen- 
tum with real spin and angular momentum of the local 
basis. This method is easily generalized to study the 
interface states between two semi-infinite crystals, e.g., 
HgTe and CdTe, which is essential in the present work 
as will be explained below. 

As is well-known for HgTe quantum wells, the con- 
finement effect along the z-direction opens a small gap 
around the Fermi level which makes the quantum well an 
insulator. By the analysis in Ref. based on the k • P 
approximation, a quantum phase transition from a topo- 
logically trivial phase to a non-trivial phase (QSH phase) 
occurs for CdTe/HgTe quantum- wells at some critical 
thickness of HgTe layer, which is signaled by a level cross- 
ing between the El and HHl subbands. We study this 
transition using the more realistic full TB Hamiltonian 
given above. First we choose periodic boundary condi- 
tions in X and y to obtain the 2d subband spectrum for 
different thicknesses Nh of HgTe layers. In Fig[l] (b), 
we plot the subband spectrum at the F-point as a func- 
tion of HgTe "half -layers. Whereas Ref. [l| predicts 
that the subbands will cross, we find an anti-crossing at 
the F-point at a critical layer thickness dc — 9a. This 
difference between the work of Ref.fl'] and these TB re- 
sults can be explained by the presence of BIA in the 
zinc-blende lattice which is ignored in the previous pa- 
per. From the k • P perspective, an additional term 
H' — Ckkz{Jz, {Jx — -^y )} is allowed in the bulk Hamilto- 
nian once the point group symmetry is reduced by BIA 
to -D2(i[ll], with Jx,Jy,Jz the spin-3/2 matrices. De- 
spite being a fcz-dependent term in the bulk system, in 
the quantum well H' generates a constant term with fi- 
nite matrix elements connecting the |i?l± > subbands 
with the \HH1^ > subbands near the F-point, which 
can be derived following the effective 2d k ■ P approach 
of Ref.[l]. This term qualitatively affects the physics in 
exactly the manner predicted by the TB model. 

Upon further investigation, we find that the inversion 
asymmetry shifts the crossing point from the F-point to 
eight non-zero A:-points around it. In Fig[T] (a), we plot 
the subband spectrum at one of these crossing points 
k = (0.017, 0.008)-^ where a clear level crossing between 
El and HHl is observed. (Notice that away from F- 
point, the inversion asymmetry removes the double de- 
generacy as expected.) The other seven crossing points 
are determined by the point-group star of this one. We 
find that even though the level crossing at F point is 
avoided, the gap-closing quantum phase transition be- 
tween two insulating regimes predicted previously in the 
simplified modelil] still exists in the more realistic case. 

After identifying the presence of the quantum phase 
transition, the next natural question is whether the topo- 
logically non-trivial phase on the thicker side of the tran- 
sition predicted by the k • P calculation will survive the 
inversion symmetry breaking. The most convincing way 
to answer this question is to calculate the edge states 
directly. By the recursion method we obtain the charge 
and spin density of states defined above. The charge den- 
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sity of states for a quantum well structure with 20 and 
10 "half -layers are plotted in FigsH] (a) and (b) respec- 
tively. Sharp peaks appear in the gap for the quantum 
well with 20 half-layers and are absent in the 10 half-layer 
system. As required by time reversal symmetry (TRS), 
the energy levels at the F point must be doubly degen- 
erate, but can be split at finite ky. The energy level 
splitting near fcj^ = as a function of ky is plotted in the 
inset of Fig. [2ja), which shows a perfect linear dispersion 
indicating a level crossing of the edge states at ky = 0. 
The spin densities of states for two fcy's with opposite 
signs are plotted in Fig. ^c) . If the chemical potential 
lies between the two peaks, only the lower branch of the 
edge states is occupied and a spin current will be car- 
ried by the edge states. With the recursion method we 
can also obtain the charge density of states on the inner 
layers away from the edge. In the inset of Fig. [2jc), we 
plot the height of the in-gap peak for fcj, = on the dif- 
ferent layers, which decays very fast from the edge and 
thus demonstrates that the in-gap peak is produced by 
the edge states. Finally, we plot the edge state disper- 
sion in Fig. IHd) with a color intensity plot generated 
from the density of states. We find that the edge states 
merge with the bulk states very quickly as ky moves away 
from the F-point. Consequently, whenever the chemical 
potential lies in the bulk gap, the only low-energy states 
crossing the Fermi level are one Kramers's pair of edge 
states. In other words, the low energy behavior of the 
quantum well in the bulk insulating region is described 
by an odd number of pairs of Id channels propagating 
on each edge. According to RefH, such a Id liquid is a 
"helical liquid," and cannot be realized in any pure Id 
system that preserves TRS. It can only exist as the edge 
theory of a 2d QSH insulator. In this way, the results 
of our calculations provide convincing evidence that the 
HgTe/CdTe quantum well is a QSH insulator character- 
ized by a non-trivial Z2 topological invariant P, [3| . 

The same recursion method can also be applied for the 
2d surface states of strained 3d HgTe, which was recently 
suggested to be topologically nontriviallT^. However, 
the TB model ([1]) applied to a semi-infinite system with a 
2d surface always generates some surface states, even for 
a trivial band insulator, e.g., CdTcflG*! Physically, these 
surface states correspond to a dangling sp^-hybridized 
bond at each surface atom, and are thus strongly depen- 
dent on the details of the surface physics, such as surface 
reconstruction and disorder. Since in the present paper 
we are only concerned about the topological properties, 
which are insensitive to the details of surface physics, we 
can choose a surface regularization that removes the triv- 
ial surface states and leaves only the topological ones. 
Therefore, we will focus on the interface between bulk 
HgTe and CdTe, where the dangling bonds of HgTe are 
coupled to CdTe so the trivial surface states vanish. Since 
CdTe can be adiabatically connected with vacuum by 
taking its band gap to infinity, the topological properties 
of HgTe/CdTe interface are determined by HgTe. 

Since bulk HgTe is a semi-metal, we need to apply a 
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FIG. 1: (a) The subband splitting as the function 
of layer thickness at one of the crossing points k = 
(0.017, 0.008)^. (Inset) Schematic diagram of all 8 crossing 
points in the 2d BZ (b)The subband splitting as the func- 
tion of layer thickness at the F point. 
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FIG. 2: (a)The density of states at the edge of the quantum 
well with layer thickness d — 20(a/2). (The inset plot shows 
the linear energy splitting of the edge states in the very small 
region near ky = 0.)(b)The density of states at the edge of 
the quantum well with layer thickness d = 10(a/2). (c) The 
spin density of states at the edge of the quantum well with 
thickness d = 20(a/2). (The inset plot shows the decay of 
density of states at the peak energy of the spectra with ky = 
0.)(d)The intensity color plot on the energy-momentum plane 
for the density of states at the edge of the quantum well with 
thickness d = 20(a/2). 
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FIG. 3: Intensity color plot in the energy-momentum plane 
for the density of states at the HgTe — CdTe interface. The 
uniaxial strain is applied along the [001] direction by choosing 
the c/a ratio to be 0.98 and the HgTe/ CdTe interface is 
chosen along [100] direction. 



small compressive strain, along say the [001] direction, to 
make it an insulator ll| . The tight binding parameters 
for the strained HgTe are obtained by fitting the LDA re- 
sults from the plane wave pseudo potential method p^. 
In the present paper, we apply the constrain along the 
[001] direction (z-direction) and locate the HgTe /CdTe 
interface along both the [100] (x-direction) and [001] di- 
rections. And for both cases, we find non-trivial surface 
states. In Fig. 3 we plot the charge density of states in 
a color intensity plot. The single pair of surface states 
is clearly seen in the bulk insulating gap and they cross 
at the F-point. There are no other surface states in the 
entire zone. Any pure 2d band theory that respects TRS 
must have the same number of Kramers's pairs on each 
of the time-reversal symmetric wave vectors because all 
of the 2N energy bands must be paired on these wavevec- 
tors. As a result, the 2d surface states of HgTe cannot 
emerge from any pure 2d surface effect, and only as a 
consequence of bulk topology. 

It is useful to try to understand this result from a con- 
tinuum k • P perspective. HgTe has a non-trivial topo- 
logical structure because the band-structure is only "in- 
verted" near the F-point. The fact that an occupied band 
at this point has F^ character means that the Z2 invariant 
picks up an extra factor of (—1) (if we ignore the small 
BIA) making it non-trivial. Due to strong orbital mixing 
the F^ character is washed out as one moves away from 
the F-point and at the other special TR invariant points 
the "inverted" structure is absent. Therefore, we should 
be able to understand the topological properties from the 
band structure only near the F-point. The key point is to 
consider the full 6-band Kane modelfl9t instead of just 
the reduced 4-band Luttinger model[20]. If we only keep 
the bands in the Luttinger model the topological struc- 
ture is absent. In the presence of uniaxial compressive 
strain along the (001) direction an insulating gap opens 
between the heavy-holc(HH) and light-hole (LH) bands 
by pushing the HH band downward in energy. For a mo- 



ment we will ignore the HH band and focus only on the 
LH and F^ (E) band. From the form of the Kane model, 
the coupling of the LH and E bands near the F-point is 
exactly a 3d anisotropic massive Dirac Hamiltonian if we 
ignore BIA and keep the leading order in k. The Dirac 
Hamiltonian preserves parity symmetry and we can la- 
bel the bands by parity eigenvalues. Since the coupling 
is linear there must be one even (doubly degenerate) and 
one odd (doubly-degenerate) band. We expect that when 
the odd parity band lies below the even band then there 
will be a non-trivial Z2 invariant which indicates an odd 
number of pairs of surface states that cross at TR invari- 
ant points [14]. The presence of the HH band will change 
the features of the spectrum but it does not change the 
presence of the surface states, or their protected crossing, 
as long as the strain induced gap is open. The system 
will remain a 3d topological insulator when the HH band 
is coupled, and when BIA terms are added, as long as the 
bulk gap does not close. To show evidence of our state- 
ments we solve the 6-band Kane model on a cylinder. 
First we solve the model with the HH band completely 
decoupled from the LH and E bands (Fig. 2^) . Here the 
HH band remains flat and is split from the LH band by 
the strain induced gap. In the gap there are clear, linearly 
dispersing surface states which traverse the gap between 
the LH and E bands. Nothing occurs at the other special 
points in the BZ, therefore this is a strong topological in- 
sulator. Turning on the coupling to the HH band changes 
features of the band structure but does not change the 
topology of the state since the gap between the LH and 
HH bands never closes. It is clear from Fig. that even 
when the HH band is fully coupled the system is still a 
strong topological insulator with surface states crossing 
at F. 
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FIG. 4: Band structure near the F-point for (a) Decoupled 
HH band (b)Full HH band coupling. Surface sates are shown 
in red. Strain induced gap is artificially large so that surface 
states are clearly visible. Generally, topologically non-trivial 
surface states exist for any finite compressive strain. 

This work is supported by the NSF under grant num- 
bers DMR-0342832 and the US Department of Energy, 
Office of Basic Energy Sciences under contract DE- AC03- 
76SF00515, the MARCO Center on Functional Engi- 
neered Nano Architectonics (FEN A), and the Knowledge 
Innovation Project of the Chinese Academy of Sciences. 



5 



[1] B. A. Bernevig, T. L. Hughes, and S.C. Zhang, Science 

314, 1757 (2006). 
[2] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 

(2005). 

[3] B.A. Bernevig and S.C. Zhang, Phys. Rev. Lett. 96, 

106802 (2006). 

[4] L. Fu, C. Kane, and E. Mele, arxiv: cond-mat/0607699. 
[5] J. E. Moore and L. Balents, arxiv: cond-mat/0607314. 
[6] C. Wu, B.A. Bernevig, and S.C. Zhang, Phys. Rev. Lett. 

96, 106401 (2006). 
[7] C. Xu and J. Moore, Phys. Rev. B 73, 045322 (2006). 
[8] S. Murakami, N. Nagaosa, and S.C. Zhang, Science 301, 

1348 (2003). 

[9] J. Sinova et. al, Phys. Rev. Lett. 92, 126603 (2004). 
[10] S. MuraJsami, N. Nagaosa, and S.C. Zhang, Phys. Rev. 

B 69, 235206 (2004). 
[11] S. Murakami, N. Nagaosa, and S.C. Zhang, Phys. Rev. 

Lett. 93, 156804 (2004). 



[12] L. Fu and C. Kane, arxiv: cond-mat/0611341. 
[13] J. Slater and G. Koster, Phys. Rev. B 94, 1498 (1954). 
[14] A. Kobayashi, O. Sankey and J. Dow, Phys. Rev. B 25, 
6367 (1982). 

[15] J. K. M. S. I. Turek, V. Drchal and P. Weinberger, Elec- 
tronic structure of disordered alloys, surfaces and inter- 
faces, Kluwer Academic publishers, 1997. 

[16] Garnett W. Bryant, Phys. Rev. B 35, 5547 (1987). 

[17] Z. Fang and K. Terakura, J. Phys.: Cond. Mat. 14, 3001 
(2002). 

[18] R. Winkler, Spin-orbit coupling effects in two- 
dimensional electron and hole systems, Springer 
Tracts in Modern Physics, 2003. 

[19] E.G. Kane, J. Phys. Chem. Solids 1, 249 (1957). 

[20] Y.R. Lin-Liu and L.J. Sham, Phys. Rev. B 32, 5561 
(1985). 



